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Abstract 

Logarithmic gaps have been used in order to find a periodic com- 
ponent of the sequence of prime numbers, hidden by a random noise 
(stochastic or chaotic). It is shown that multiplicative nature of the 
noise is the main reason for the successful application of the logarith- 
mic gaps transforming the multiplicative noise into an additive one. 
A relation of this phenomenon to spontaneous neuron activity and to 
chaotic brain computations has been discussed. 

1 Introduction 

The prime number distribution is apparently random. The apparent ran- 
domness can be stochastic or chaotic (deterministic). It is well known that 
the Riemann zeros obey the chaotic GUE statistics (if the primes are inter- 
preted as the classical periodic orbits of a chaotic system, see for a recent 
review [1]), whereas the primes themselves are believed to be stochastically 
distributed (Poissonian-like etc. [2]). However, recent investigations suggest 
that primes themselves "...could be eigenvalues of a quantum system whose 
classical counterpart is chaotic at low energies but increasingly regular at 
higher energies." [3]. Therefore, the problem of chaotic (deterministic) be- 
havior of the moderate and small prime numbers is still open. Moreover, there 
are also an indication of periodic patterns in the prime numbers distribution 
[3]-[7]. These patterns, however, has been observed in probability distribu- 
tion of the gaps between neighboring primes and not in the prime numbers 
sequence itself (see also Ref. [8]). The intrinsic randomness (stochastic or 
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chaotic) of the primes distribution makes the problem of finding the peri- 
odic patterns in the prime numbers themselves (if they exist after all) a very 
difficult one. 

A physicists may ask: Why should one be interested in finding these 
patterns? The answer is: comparison. If one can recognize patterns in 
an apparently random system, then one can compare these patterns with 
the patterns known for some other systems of interest. We have already 
mentioned the comparison with certain quantum systems. Another intrinsic 
comparison can be made with the computational properties of brains, where 
the natural numbers certainly should play a crucial role. The neuron signals 
are also apparently random. Can one compare patterns observed in these 
signals with the patterns in the prime numbers sequence in a constructive 
way, in order to shed a light on the computational apparatus of the brains? 

2 Logarithmic gaps 

It is believed that properties of the gaps between consecutive primes can 
provide a lot of information about the primes distribution in the natural 
sequence. The so-called prime number theorem states that the "average 
length" of the gap between a prime p and the next prime number is propor- 
tional (asymptotically) to Inp (see, for instance, Ref. [2]). This implies, in 
particular, statistical non-stationarity of the prime numbers sequence. That 
is a serious obstacle for practical applications of the statistical methods to 
this sequence. For relatively large prime numbers one can try to overcome 
this obstacle by using relatively short intervals [2 J . Another (additional) way 
to overcome this obstacle is to use logarithmic gaps (logarithms of the gaps). 
The non-stationarity in the sequence of the logarithmic gaps is considerably 
'slower' than that in the original sequence of the gaps themselves. The log- 
arithmic gaps have also another crucial advantage. If the noise corrupting 
the gaps sequence has a multiplicative nature, then for the logarithmic gaps 
this noise will be transformed into an additive one. It is well known that 
the additive noises can be much readily separated from the signal then the 
multiplicative ones (see below). 

Before starting the analysis let us recall two rather trivial properties of 
the gaps, which will be used below. The sequence can be restored from the 
gaps by taking cumulative sum. Pure periodicity in a sequence corresponds 
to a constant value of the gaps (the period). 
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Figure 1: Autocorrelation of the v(n) functions for the ln-sequence (circles) 
in the interval 9 x 10 4 < n < 10 5 . The solid line correspond to autocorre- 
lation function computed for the telegraph signal generated by the Rossler 
attractor fluctuations overcoming the threshold x — 7. In order to make 
the autocorrelation functions comparable a rescaling has been made for the 
Rossler autocorrelation function. 

Let us take, as a first step, logarithms of the gaps for the sequence of 
prime numbers. Then, let us compute cumulative sum of these logarithms as 
a second step. Then, we will multiply each value in the cumulative sum by 
10 and will replace each of the obtained values by a natural number which 
is nearest to it. As a result we will obtain a sequence of natural numbers: 
7,14,28,35,49,55,.... which will be called as ln-sequence. 

Let us now define a binary function v(n) of natural numbers n — 2, 3, 4, 
which takes two values +1 or -1 and changes its sign passing any number 
from the ln-sequence. This function contains full information about the ln- 
sequence numbers distribution in the sequence of natural numbers. 

Let us consider a relatively short interval: 9 x 10 4 < n < 10 5 . For this 
interval the autocorrelation function 

C(n, r) = (v(n)v(n + r)> - (v(n))(v(n + r)> (2.1) 

computed for the v(n) function of the ln-sequence will be approximately 
independent on n. Figure 1 shows the autocorrelation function (circles) 
computed for the v(n) function of the ln-sequence for the 'short interval': 
9 x 10 4 < n < 10 5 . Figure 2 shows an autocorrelation function (circles) com- 
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Figure 2: As in Fig. 1 but for interval 1000 < n < 2 x 10 4 . 

puted for the v(n) function of the ln-sequence for interval of much smaller 
natural numbers: 1000 < n < 2 x 10 4 (here we rely on the much slower non- 
stationarity in the logarithmic gaps sequence in comparison with the original 
gaps sequence). 

3 Models 

In order to understand origin of the oscillating autocorrelation function shown 
in Figs. 1 and 2 let us consider a very simple telegraph signal, which allows 
analytic calculation of its autocorrelation function. The telegraph signal 
takes two values +1 or -1 and it changes its sign at discrete moments: 

t n = nT + ( (3.1) 

where ( is an uniformly distributed over the interval [0, T] random variable, 
n = 1,2,3... and T is a fixed period. If q is a probability of a sign change 
at a current moment (0 < q < 1), then the autocorrelation function of such 
telegraph signal is: 

C(r) = (n - r/T)(2q - l)"" 1 + (r/T - (n - l))(2q - If (3.2) 

in the interval (n — 1)T < r < nT. Figure 3 shows the autocorrelation 
function Eq. (3.2) calculated for q = 0.25, as an example (cf. with Figs. 1 
and 2). In the model Eq. (3.1) ( was taken as a pure stochastic variable. 
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Figure 3: Autocorrelation function for the simple model telegraph signal: 
Eq. (3.2) with q = 0.25. 
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Figure 4: X-component fluctuations of a chaotic solution of the Rossler sys- 
tem Eq. (3.3) (a = 0.15, b = 0.20, c = 10.0). 
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This variable, however, can be also a chaotic (deterministic) one. Let us 
consider a chaotic solution of the Rossler system [5] 



dx 
~dt 



{y + z); 




dz 

— = b + xz — cz 



(3.3) 



where a, b and c are parameters (a reason for this system relevance will be 
given in Section 4). Figure 4 shows the x-component fluctuations of a chaotic 
solution of the Rossler system. 

Let us consider a telegraph signal generated by the variable x crossing 
certain threshold from below (Fig. 4). This telegraph signal takes values: 
+ 1 or -1, and changes its sign at the threshold-crossing points. Then, the 
fundamental period of the Rossler chaotic attractor provides the period T in 
Eq. (3.1), whereas the chaotic fluctuations of the variable x of the attractor 
provide chaotic variable ( in the Eq. (3.1) for the considered telegraph signal. 
Autocorrelation function computed for this telegraph signal has been shown 
in Figs. 1 and 2 as the solid line. In order to make the autocorrelation 
functions comparable a rescaling has been made for the Rossler telegraph 
signal's autocorrelation function: the scaling coefficient is equal to 0.36 for 
Fig. 1 and to 0.30 for Fig. 2. 

4 Prime numbers and spontaneous neuron ac- 
tivity 

All types of information, which is received by sensory system, are encoded 
by nerve cells into sequences of pulses of similar shape (spikes) before they 
are transmitted to the brain. Brain neurons use such sequences as main in- 
strument for intercells connection. The information is reflected in the time 
intervals between successive firings (interspike intervals of the action potential 
train, see Fig. 5b). There need be no loss of information in principle when 
converting from dynamical amplitude information to spike trains and the 
irregular spike sequences are the foundation of neural information process- 
ing. Although understanding of the origin of interspike intervals irregularity 
has important implications for elucidating the temporal components of the 
neuronal code the problem is still very far from its solution. 

In a recent paper [TU] a membrane threshold-crossing model of sponta- 
neous (without external stimulus) neuron firing based on the chaotic Rossler 
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Figure 5: Mapping of a spike train (figure lb) into a telegraph signal (figure 
la). 
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Figure 6: Autocorrelation functions for the telegraph signals corresponding 
to a spontaneous activity of a hippocampal (CA3) singular neuron (solid 
curve) and to the spike train generated by the Rossler attractor fluctuations 
overcoming a threshold (circles). In order to make the autocorrelation func- 
tions comparable a rescaling has been made for the Rossler autocorrelation 
function. The data has been taken from the Ref. |10|. 
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Figure 7: Autocorrelation of the v(n) functions for the ln-sequence (circles) in 
the interval 9x 10 4 < n < 10 5 . The solid curve corresponds to autocorrelation 
function computed for the telegraph signals corresponding to a spontaneous 
activity of a hippocampal (CA3) singular neuron (cf. Fig. 3). In order to 
make the autocorrelation functions comparable a rescaling has been made 
for autocorrelation of the v(n) function. 

attractor was suggested. Good agreement between the model and the spik- 
ing time-series, obtained in vitro from a spontaneous activity of hippocampal 
(CA3) singular neurons (rats brain slice culture), was reported in this paper. 
The spiking time-series were mapped into a telegraph signal as it is shown 
in Fig. 5. The comparison between the experimental data and the signal 
generated by the Rossler attractor fluctuations overcoming a threshold (cf. 
previous Section) is shown in Fig. 6. 

The hippocampus is a significant part of a brain system responsible, in 
particular, for spatial memory, body orientation and navigation. These func- 
tions of brain are related to the very sophisticated real-time computations, 
which the brains make with astonishing effectiveness. One should distinguish 
between the inner brain computations and our conscious arithmetic calcula- 
tions, which are performed in the more specialized brain's departments such 
as the intraparietal sulcus [11]. Therefore, if one wants to investigate how 
the inner brain computational abilities are related to neuron activity the hip- 
pocampal pyramidal cells are an appropriate subject for this investigation. 
Spontaneous activity in brain slice preparations purely reflects the intrin- 
sic properties of local circuits and individual neurons. The comparison of 
the spontaneous activity of the hippocampal singular neurons with the ran- 
dom properties of the sequence of prime numbers can shed certain light on 
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Figure 8: As in Fig. 7 but for interval 1000 < n < 2 x 10 4 . 

the long standing problem: whether mathematics belongs completely to the 
outer world or it is a production of our brain activity (at least in the form 
of the natural numbers). For instance, in order to work together the brain 
neurons have to make an adjustment of their rhythms. The main problem for 
this adjustment is the very noisy environment of the brain neurons [T2]. For 
pure periodic inner clocks this adjustment would be impossible due to the 
noise. Nature, however, has another option. This option is a chaotic clock 
(see Ref. [10] and references therein). In chaotic attractors certain character- 
istic frequencies can be embedded by broad-band spectra, that makes them 
much more stable to the noise perturbations. The Riemann's hypothesis, 
for instance, which states that the apparent randomness in the distribution 
of prime numbers only comes from a noise (otherwise the prime numbers 
are as regularly distributed as possible), can reflect an optimization of such 
mechanism. 

On the other hand, the compatibility of the statistical properties of the 
neuron spontaneous activity with the statistics of the prime numbers (see 
Figs. 7 and 8) could be a key for understanding of the inner computational 
methods of the brains. If brains use chaotic dynamics (Fig. 6) also for 
computational processing [13], then such compatibility should play a crucial 
role in the inner brain computations. The reduction to the additive random 
processes by means of the logarithmic gaps (see previous Section) could be 
a part of these inner brain computations. If the gaps can be considered as 
a result of a multiplicative random process, i.e. each gap between consecu- 
tive primes can be represented as a product of independent random variables 
(the random variables could take their values from the prime numbers set, 
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including 1), then logarithmic gap represents the sum of random variables 
(i.e. an additive random process). It is known that unlike the additive ran- 
dom processes the multiplicative random processes are much more complex 
and difficult for analysis. 

I thank Y. Ikegaya for sharing the neuron data and discussions. 
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